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We present a new semiclassical technique which relies on 
replacing complicated classical manifold structure with sim- 
pler manifolds, which are then evaluated by the usual semi- 
classical rules. Under circumstances where the original mani- 
fold structure gives poor or useless results semiclassically the 
replacement manifolds can yield remarkable accuracy. We 
give several working examples to illustrate the theory pre- 
sented here. 



I. INTRODUCTION 

Semiclassical methods are based on classical mechan- 
ics; the relevant classical manifolds form the "skeleton" 
to which the wave function is attached. However, associ- 
ating a semiclassical wave function with a classical man- 
ifold is problematic when enclosed loops in phase space 
have area less than h. If that happens, typically two or 
more stationary-phase points corresponding to distinct 
contributions to the semiclassical wave function are not 
well separated in phase. Consequently, the stationary- 
phase (SP) method, on which semiclassical approxima- 
tion is based, breaks down. Worse, in some important 
situations, such as the universal homoclinic oscillations 
associated with chaotic regions of phase space, the bad 
enclosed loops may be repeated many times in a small 
region 0. 

Although the semiclassical wave function diverges near 
caustics lurking near the small area loops, the exact 
quantum wave function is of course well-behaved. This 
is one manifestation of quantum smoothing over classi- 
cal detail. We shall quantify below how this smoothing 
occurs. We suggest and test a new way of looking at the 
smoothing process in terms of "replacement manifolds" , 
in which new, well behaved classical manifolds are sub- 
stituted for the original badly behaved ones. One ap- 
plication of this approach is uniformizing a semiclassical 
wave function in the vicinity of the ubiquitous homoclinic 
oscillations of a chaotic system. 

It is widely understood that very small changes in po- 
tentials or walls of billiards can have little effect on quan- 
tum eigenstates, but can drastically affect the classical 
manifold structure. This immediately implies a many- 
to-one relationship between the classical manifolds and 
quantum wave functions. That is, many different under- 
lying classical manifold patterns correspond to the same 
wave function. Now, whatever the choice of manifold 
(amongst these equivalent forms) we can suppose an ap- 
propriate uniformization exists which gives the correct 



wave function. However, the non-uniformizcd, simple 
semiclassical limit for the various choices of manifolds 
can differ drastically in accuracy. It is this fact which we 
exploit in the present paper. 



II. STANDARD UNIFORMIZATION: THE AIRY 
FUNCTION 

The paradigm example of uniformization is the Airy 
function, i.e. the energy eigenfunction of the linear ramp 
potential. We take the potential to be of the form 



V(x) = -fix 
The Schrodinger equation is 

h 2 •■ 
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il> E {x) - (px + E)^ E (x) = 



After the shift 

1>e(x) = Vo(a: + £//?) 
and scaling transformation 

_ to 1 / 3 /? 1 / 3 

U ~ fi2/3 X ~ 1 X 

the Schrodinger equation reads 

+ v$>(u) = 0; 
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with <3>(w) = t/Jo(x). Eq. (||) is exactly solved by going to 
the momentum representation. We have 



l> 2 -%)*"" " 



with solution 
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Thus 



'$(p) dp 



cos(p 3 /6 — p u) dp 



>i/2 Ai f- 2 ^u\ 



(8) 



1 



and 



if)(x) = T) Ai(-2 1/3 7 x) = 77 Ai ^> 
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SP approximation for the momentum integral gives 
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Noting that 
we get finally 



(2m) 3 / 2 



(9) 



(10) 



: COS 



J U p(u')du'-^ 



This is the standard WKB form for the linear ramp po- 
tential. Figure [l] illustrates the situation for the semi- 



Area ~ classical /fvj 
probability 




FIG. 1. Diagram showing the areas in phase space impor- 
tant to the semiclassical approximation to the Airy function. 

classical position wave function, which becomes inaccu- 
rate when the area enclosed by the vertical a>line and the 
energy contour corresponding to the state \E) falls below 
Planck's constant. Upon entering the "bad" region, the 
original momentum integral may be substituted for the 
WKB result, giving a uniformization. We have an in- 
tegral over "initial momentum," making this example a 
type of "initial value representation" (IVR) as first dis- 
cussed by Miller |^,|| . Of course in this case we know the 
result of the integral is an Airy function, but in general 
there may be integrals which are not known analytically 
or even close to ones that are known. 

In this paper we encounter more complex situations 
which arc induced by nonlinear interactions, including 



chaos, reflection from corrugated surfaces, etc. We will 
show that intricate manifolds can sometimes be replaced 
by much simpler ones, which may themselves be evalu- 
ated semiclassically in any basis, using a much simplified 
uniformization if required. 



III. REPLACEMENT MANIFOLDS 

We now introduce an example of such a problematic 
Lagrangian manifold, defined by a position-momentum 
relation 



P(x) = po 



■ hXax sin ax 



(11) 



This model was first introduced in Ref. [Q but was not 
explored there. The function ( pd| ) has oscillations of in- 
creasing amplitude and frequency, which however have 
the same area J p(x)dx = Xh between successive zeros 
(see Fig. |^, small value of A corresponds to small area 
of loops) . This is also a characteristic property of homo- 
clinic oscillations near a periodic orbit. The momentum 
p(x) spans an ever larger range as x increases, but we 
now show that for A < 1 the manifold's semiclassical 
behavior can be understood by replacing it with three 
smooth manifolds. The rules for replacement are simple 
to derive. 

The action function S(x) is the integral of p(x): 

f x 1 
S(x) = / p(x')dx' =p x + -h\(cosax 3 - l) . (12) 
Jo 2 



The "wave function" 

iP(x)=A(x) e tS{x)/h 
can be approximated by 

i/j(x) m A(x) e - 4A /2 e *P<Wft ^i + cos(ax 3 ) 
if>(x) m A{x)e- lX/2 e tpoX / h 



(13) 



(14) 



(15) 
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if A < 1. For small A this expression is nearly iden- 
tical with ([[3]) with the original action (p^), because 
cos(ax 3 ) is bounded by ±1. However semiclassically it 
has the interpretation of the sum of three smooth clas- 
sical manifolds, i.e. p — po, p = po + 3hax 2 and p — 



Po — 3 hax 2 . They have the weights A(x), i\e lX l 2 A{x), 
ije~' tX / 2 A(x), respectively. The situation is depicted in 
Fig.|. 

Almost every discussion of the relation of classical and 
quantum mechanics for chaotic systems alludes to quan- 
tum smoothing, but here we have seen explicitly one way 
this smoothing comes about. For all reasonable purposes 



2 



the three smooth classical manifolds accurately replace 
the rapid oscillations of the original manifold. Note too 
that depending on the parameters the outlying manifolds 
p = po + 3 hax 2 , and p = po — 3 hax 2 lie far beyond 
the limits of the original distributions. We shall discuss 
interesting consequences of this in the next section. 

So far, we have considered the wave function only in 
position representation where caustics are absent and the 
semiclassical form ( |l3| ) was accurate. The situation dete- 
riorates drastically in momentum space. In order to find 
the wave function in momentum representation, we have 
to sum over all contributions from intersections of a hor- 
izontal line (corresponding to a momentum eigenstate) 
with the oscillating manifold ([ll]). When A < 1, the ad- 
jacent intersections will be separated by a phase AS / h 
that is smaller than 1 for any classically allowed momen- 
tum, and therefore the standard semiclassical approxi- 
mation will break down for all classical momenta). Usual 
Airy-type uniformization methods of Section || (also see 
HD) do not serve under these circumstances, because 
the manifold is dense with caustics (see the filled-in ar- 
eas A Si and AS2 in Fig. ||). Apparently, we have no 
other choice but to part with semiclassical method and 
use a numerical Fourier transform (IVR) over the po- 
sition wave function to obtain 4>(p). However using the 
replacement manifolds (RMs) is a much simpler and more 
intuitive approach which can rescue the situation. To find 
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FIG. 2. Areas in phase space important to mo- 
mentum wave function associated with the manifold 
p(x) = po — ^haXx 2 sin(ax 3 ), X — 0.5. Primitive semiclas- 
sical approximation fails since both areas A Si and AS2 can 
be smaller than ft. 

the momentum wave function cj>{p) associated with the 
original manifold we add (with appropriate weights) wave 
functions 4> n {p) corresponding to the three RMs. Each of 
these partial wave functions is equal to a Fourier trans- 
form of a corresponding position wave function ip n (x), 

MP) = ~i= I dxi, n {x)e-^' h . (16) 

Except for a small region close to po (highlighted in 
Fig. ||), this integral can be evaluated by SP method! 
In other words, although the semiclassical approximation 
completely fails for the original manifold, it works almost 
everywhere for the replacement manifolds. Moreover, as 
we show in the following section, the standard Airy uni- 
formization procedure of Section O can be exploited to 



correct the inaccuracy in a narrow region near pq. We will 
see that the momentum wave function found by applying 
SP approximation to RMs is in excellent agreement with 
numerical solution. We emphasize that the RMs do not 
yield the original, badly behaved semiclassical result, but 
rather something much more accurate and much simpler. 




x 



FIG. 3. Original (solid line) and replacement manifolds 
(dashed lines) for p(x) = po — fftaAa; 2 sin(aa; 3 ), A = 0.5. 



IV. UNIFORMIZATION WITH REPLACEMENT 
MANIFOLDS IN A GENERALIZED MODEL 

In order to capture other properties of homoclinic oscil- 
lations in our model, we consider a generalized manifold 

p(x)= Po -h\^l S max), (17) 

in which ax 3 from ( pT| ) is replaced by a generic smooth 
function £(x). The presence of the derivative en- 
sures equal loop areas between successive zeros of £ (x) . 
Example of a manifold with £(x) — alog^ is displayed 
111 Fig. |. Manifold |TJ) can be obtained from a hori- 
zontal manifold representing a momentum state |po) by 
a canonical transformation of coordinates |jj 

1 

x = Xq, p = po - -%\- — sin£(a;o), (18) 
2 axo 

generated by a function 

F 2 (x,p ) =po%+~ ftX cos £(x), (19) 

We find the uniform semiclassical transformation ele- 
ment (p\po), and thereby fix the amplitude prefactor to 
be A(x) = {2ttK)^ 1 / 2 ; this however retains all the prob- 
lematic features of the manifold. We solve the problem 
to all orders in A by replacing the original manifold with 
an infinite series of RMs. First we note that the semi- 
classical transformation element 

(x\p ) sc = -L=e^^' h (20) 
V27m 

is fairly accurate since there are no caustics in this repre- 
sentation (to each final x corresponds a single initial xq). 
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The momentum-representation element is obtained by a 
Fourier transform, 



(p|po)~-= / dxe-^/ h (x\ Po ) sc 
V 27m j 



(21) 



If evaluated by SP method, the integral yields a very 
poor semiclassical result since there exist many coalesc- 
ing SP points. An accurate uniform answer is obtained by 
evaluating integral ( f2l| ) exactly. The accuracy is usually 
further improved by changing the integration variable 
from x to #0 (in our case both forms are equivalent since 
x = Xq). This uniform version of (p\po) is in IVR form 
|§,|J, but can best be evaluated by writing (x\po) as a 
sum over RMs. Recognizing that the factor e l ? cos t( x ) in 
(x\po) sc = (27rfi) -1 / 2 exp [ip x/H + cos£(a;)] is a gen- 
crating function for Bessel functions, we can extend the 
sum from previous section beyond the first order in A. In 
fact we obtain an infinite sum convergent for any A, 



(x\Po) 



"I" t A (5) 




jn^ntix) (22) 



where (x\po) n = (2nh) 1 / 2 exp {i [pox + nH£(x)] /%}. We 
can rewrite ( pl| ) as 

(p|po)= £ J n [ dxe-^ h (x\p Q ) n (23) 



A. The replacement-manifold method is therefore not re- 
stricted to the regime where loop areas are smaller than 
h. It also works in the strongly chaotic regime where 
A > 1 (and where the standard semiclassical approxima- 
tion holds) if all manifolds up to \n\ ~ A arc included 
in the sum. This follows since J n (A/2) considered as a 
function of n decays exponentially fast for \n\ > A/2. In 
other words, we only need to include several RMs beyond 
those intersecting the original manifold. Physically, the 
manifolds outside of the range of the original one lie in 
a classically forbidden region in which the wave function 
is exponentially suppressed. 

In opposite case (A < 1) which interests us the most, 
considering RMs up to |n| = 1 will suffice while the 
"primitive" semiclassical approximation fails completely 
even for p — pq. In this regime A can be thought of as a 
parameter defining the strength of a perturbation which 
causes the manifold to oscillate around a simple manifold 
that describes the system if perturbation is absent. 

To complete our solution for a specific function £(x), 
we must evaluate (p\po) n - For n — 0, we cannot use SP 
approximation (because the action is linear), but inte- 
gral (p3) is trivial and we obtain in general 



(Pbo)o = 5(p-Po) 



(26) 



For n 7^ 0, SP method applied to (gj) yields a result in 
the form of a sum over SP points x sp , 



(p\Po)n. 

x E 





1 




c 






d 2 fn 


-1/2 




dx 2 


exp 


K 



(27) 



in 4- — 

h 



d 2 f n 



4 Sgn dx- 



n— — oo 



E J n(~)* n (p\P° 



where 

(p\Po)n = — ^ J dx exp j I [(po - p)x + nh^(x)} j (24) 

may be identified semiclassically as a transformation el- 
ement corresponding to the n-th RM 



Pn(x) =Po + nh 



dx 



(25) 



This manifold for most functions £(x) of interest contains 
no caustics, and consequently (p\po)n allows evaluation 
by SP method! As promised, we have expressed the uni- 
form version of (p\po) as a weighted sum over semiclassi- 
cal replacement manifolds, 



(p\po, 



unif 



E 



Jn 



2 ) '"G'h'>».. 



If applied to a well-behaved (i.e. decaying fast enough at 
±00) function of po, the resulting sum converges for any 



for 



3nhax 2 has two SP points 



where f n {x,p,po) = (po ~ p)x + nh£_{x) and ^ = 



For manifold (|1 1|) , each RM p n (x) = p a 



sp 



± 



P-Po 



1/2 



3nha 

whose contributions add to give 

(p\po)n,sc = *- 1/2 [3 nh 3 a(p - p )}- 1/4 

P-Po 



(28) 



(29) 



(na)V2 



3/2 



RMs still contain one caustic at p — po, which can be 
easily uniformized by evaluating integral (p3) exactly M , 



(p\P0)n 



unif 



(3|n|a)V3ft 



Ai 



(sgn n)(p-p ) 
(3|n|a)V3fi 



Having found (p\po)n, we can calculate (p\po) from ( p3] ) 
to all orders in A. (p\p )i.sc and (p\po)i,unif agree very 
well except for a small region near po (similarly as in 
Fig. Q). Fig. H demonstrates the excellent agreement 
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Im<plR)> 




P-Po 



FIG. 4. Momentum wave function: comparison of numer- 
ical IVR (solid line) with the RM sum up to |n| = 1 (0(A), 



dashed line) for original manifold with £(x) 
and a = 1. 



ax 3 , A = 0.5 



of RM method with the direct numerical computation of 
(p\Po) using fast Fourier transform (IVR). 

Here the replacement-manifold approach should be 
compared to the semiclassical perturbation (SCP) ap- 
proximation of Miller and Smith |^,|l^] who used pertur- 
bative classical dynamics to calculate the semiclassical S 
matrix. They first applied this method to the collinear 
collision of an atom with a diatomic molecule The 
starting point is the initial value representation of the 
classical scattering matrix, 



h 1 lUb 



9qf(qi,rii) 



1/2 



(30) 



x exp(z{$(q l ,n i ) + q f (qi,n t ) n*) - n f }}) 

where h = 1, rij, nj , qi, qf denote, respectively, the initial 
and final values of the action and angle variables describ- 
ing the internal degree of freedom, and $ is the action 
integral. For details, see Ref. [^) and ||. If the classical 
quantities listed above are calculated using the first-order 
perturbation dynamics, the S matrix takes the form 



2tt 



dq l exp {-i [(rif - m) qi + A(qi, rii)]} 



where $o 1S twice the phase-shift for the unperturbed 
potential and A{qi 1 rii) is the time integral of the per- 
turbation potential along the unperturbed trajectory Q . 
For various physical systems (such as ion-dipole collisions 
H or atom-surface scattering jL0|), A(qi,rii) has a sinu- 
soidal dependence on qi as long as perturbation remains 
small. In this regime we expect the SCP and RM approx- 
imations to give the same answer. We shall verify that 



for the scattering by a corrugated wall in Section VII 



The SCP approximation of Miller and Smith presently 
appears to be more widely applicable since it does not 
require any special property of the manifold. On the 
other hand, advantage of the RM method lies in the fact 
that it is not subjected to the smallness of perturbation. 
Provided that the full action falls into one of the classes 
discussed above, we can apply the RM method without 
approximating classical dynamics. 



V. MODEL OF HOMOCLINIC OSCILLATIONS 

Homoclinic oscillations in chaotic systems have an- 
other characteristic property: the amplitude of oscilla- 
tions (i.e. "height" of loops) increases exponentially as 



we approach an unstable periodic orbit (see Fig. ||). To 
P-Po 




FIG. 5. Original (solid line) and replacement manifolds 
(dashed lines) for £(x) = a log ^, A = 0.5. 

keep area of loops constant, the width of loops has to de- 
crease accordingly. Because systems in nature are usually 
bounded, the exponentially growing loops must eventu- 
ally twist in a complicated way to fit into available phase 
space. For sake of simplicity, we consider a model that 
satisfies the requirement of exponential growth, but is 
unbounded. Such a model can be found among general- 
ized manifolds of previous section. Specifically, we want 
to find a manifold whose homoclinic points (i.e. inter- 
sections with line p — po in our case) form a geometric 
series, 



(31) 



where A+ is the Lyapunov exponent. In other words, 
we are looking for a function £(x) such that t;(x n ) = 
2nn. The logarithmic function comes to mind first since 
log ^ = n\ + . With the correct prefactor a = 
find the desired function 



we 



£(x) = alog- 



(32) 



because £(x„) = alog^ = j^n\ + — 2nn. Substituting 

£(x) from ( p2| ) into the general form (p"7|), we obtain a 
manifold described by 



1 h\a . x . 

p = p a - cos(alog-) 

2 x a 



(33) 



and displayed in Fig. |j| Replacement manifolds (19) be- 
come 

nha 



p n {x) =p + 
In the present case, function /„ from (p?|) becomes 



fn(x,p,p ) = (po~ p)x + nha\og-, 

a 



(34) 



and its single SP point, 



nha 
P-Po' 



(35) 



For p > pq, only manifolds with n > contribute, and 
for p < Po only n < is allowed. Defining j = \n\, for 
p ^ po we replace the sum over n by sum over all positive 



5 



j. In both cases, the j-th manifold gives a semiclassical 
contribution 



(p\P< 



0/j,sc 



2tt 



1/2 



1 



\p~Po\ 
x exp 1 i sgn(p - p ) 



(36) 



j a log 



jha 



ae\p-p a \ 



with j — \n\ > 0. In this case caustics are missing be- 
cause there is a single SP point x sp . We have a finite 
integration limit at zero, but it is separated from x sp by 
an infinite action (equal to the area delimited by lines 
x = 0, p = po, and the ra-th RM, see Fig. ||). In fact, it 
can be shown that further terms in asymptotic expansion 
of (p\po)j,sc have the same dependence on fi and p — po as 
(pq), and only differ in their dependence on ja. Luckily, 
these claims can be easily verified by evaluating integral 
( pif ) analytically. This oscillatory integral is made con- 
vergent by displacing momentum p with an imaginary 
infinitesimal term — ie, e > 0, 



(p\po)n,unif = 7TT 

Znn 



(37) 



dx exp 



— (p — ie — po)x + nhalog ■ 



In this form, the answer is found by rotating the contour 
about the origin of complex plane by — -| sgn(p — p ). Let 
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FIG. 6. Contour for evaluation of {p\po)n,unif in the case 
£(a;) = a log x/a. 

us explicitly solve the case p > pq. Fig. || shows a contour 
in complex plane enclosing a region with no singularities, 
and as a result the sum of integrals along appropriate 
parts of the contour is zero, 1\ + 7r + 12 + I r = 0. In the 
limit R — > 00, r — > 0, 1% becomes our desired integral j3^ ) 
and both Ir and I r vanish, implying I\ = —Ii. Contour 
for integral I2 is the set of complex points z = — ix — 
e~'™l 2 x, x > 0. Consequently, 



(38) 



-h = ^- l2 x 
27m 

10 ( i xe~ 

dx exp < — [— (p — po — ie)(—ix) + nha log 

h 



Transforming to a dimensionless variable y = ( p Po ^ iSig j 
the integral becomes 



h = ^-r exp |(na - i)— - ma log [{p - po — ie)a/h]j 
dye~ y y ma . 



p-po-ie Jo 

Recognizing the remaining integral as a gamma function 
of complex argument 1 + ina, 



-h = 



r(i 



2tt(p - po - ie) 
x exp |(na — i)— — ina log [{p — po — ze)a/fi,]| . (39) 

For p < po, we need to close the contour in the upper half- 
plane, but otherwise the procedure is analogous. Com- 
bining the two cases, we obtain a single analytic expres- 
sion for the uniform transformation element, 



(p\po) 



xexp 



r(l + ina) 
2vr(|p — po| - ise) 

1" . - (\p - po\ - ise)a 



(40) 



(not — *)^ s — ma log 



where s = sgn(p — po). As promised, the uniform and 
semiclassical forms have identical dependence on h and 
p — Po- The only difference lies in the dependence on 
na, in particular the semiclassical form is obtained if we 
keep only the first term in the asymptotic expansion of 
r(l + ina) in (p\po) n ,unif ■ However, for our conservative 
choice of a w 33 the agreement between (p\po) n ,sc and 
(p\po)n,unif is such that they may be used interchange- 
ably for any practical purposes. Finally, in Fig. @, (p\p ) 
is applied to an initial Gaussian wave packet centered 
around Xq, Po, namely 



^*(po) 



nh 2 



-1/4 
exp 



-{P -p )Xo 



a 2 (p -Po) 2 



2h 2 



The final position wave function is 



exp 



j:F 2 (x,Po) 



and the uniform momentum wave function, 



6 / 

u n 1 



f (p) = Jo 
A 



5 )■■>>) 



* X (p\Po 



(27m) 1 / 2 

(^2)1/4 eX P 



(41) 



2a 2 



with Xgpj = jha/\p-p \ and (p|po)i,sc given by (||). 
The large oscillation near p = po in Fig. [7| corresponds 
to the zeroth RM contribution (fi l (p), the smaller wavelet 
to the right corresponds to the remaining RMs. The 
figure confirms the excellent accord between RM method 
and numerical IVR evaluation by FFT pj| . We have 
thus succeeded in uniformizing something with many of 
the properties of homoclinic oscillations near an unstable 
periodic orbit using replacement-manifold "technology" . 
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FIG. 7. Momentum wave function: comparison of numer- 
ical IVR (solid line) with the RM sum up to j = 1 (O(A), 
dashed line) for original manifold with £(x) = alog(^), 
A = 0.5 and a m 33. 



VI. COMPARISON WITH THE 
STATIONARY-PHASE METHOD 



TABLE I. Comparison of numerical, RM, and SP evalua- 
tion of |</Hp)| 2 f° r £,( x ) ~ ax > A = 10. 



n 


1 2 3 4 5 


\a n \ (num.) 


0.1776 0.3276 0.0466 0.3648 0.3912 0.2611 


K| (RM) 


0.1776 0.3276 0.0466 0.3648 0.3912 0.2611 


K| (sp) 


0.1704 0.3324 0.0343 0.3622 0.4312 oo 



S(x,p,p ) = -px + F 2 (x,p ) = (p -p) x + -HXcosax. 



Using 



15[ ) and defining r = p h ^° , we get 



1 oo 

(p\p )sc = j-F(X,r) J2 e 2 



(46) 



It is instructive to check explicitly if the expansion in 
terms of RMs agrees with semiclassical method applied 
to the original manifold when areas of loops are larger 
than h (A > 1). Let us choose an analytically solvable 
example with = ax. Corresponding manifold 



p{x) = Pa 



1 



h\a sin ax 



(42) 



has another advantage compared to manifolds in Fig- 
ures H and M. Unlike those, for A 3> 1 and small enough 
p — po, manifold ( [42| ) has all caustics in a safe distance. 
The RMs are horizontal lines p n = po+nha, independent 
of x and corresponding transformation elements can be 
evaluated exactly as (p\po) n — $(p — Po — nha). Using 
(p3|), the uniform expression for (p\po) is 



(p\po) 



unif 



oo 

E 

n= — oo 



i n 5(p-pa-nha), (43) 



so the wave function is determined by its Fourier coeffi- 
cients c„(A) = J„(A/2) i n , which is natural since manifold 
(fl2|) is exactly periodic. Moreover, each RM contributes 
only to a single momentum {p n ). Put differently, we have 
calculated the wave function for any specific momentum 
to all orders in A. Now let us find the semiclassical form 
of (p\po)- For a given momentum p there is an infinite 
number of SP points which occur in pairs (x n ,y n ), 



2tt 

x H n, y n 



TV 

a 



2n 



x 



-n, 



(44) 



where xo — — ^ arcsin xh ^ ■ Analogously to fl27h the 
semiclassical transformation element for manifold (4J) is 



(p\Po)sc = 



1 



E 



n— — oo 



d 2 s 



dx 2 



-1/2 



(45) 



x < exp 



- h S{x n ) - i- 



'X — X n 

7T 



exp 



with 



where 
F(A,r)= 



7T 



t/2 



(47) 



: )+ (f) 



1/2 



[(|) 2 -^] 1/4 

Employing the Poisson summation formula 

oo oo 

e 27rmr = S ( m - r ) 

— oo m— — oo 

and reverting to p = pq + rha, we obtain 

1 oo 

(p\po) sc = j-F(\,r) § ( n ~ r ) 

Tl— ~~ OO 

1 

= E F(Kn)6(n-r) 

n— — oo 

oo 

= Y F(\,n)6(P -p Q -nha). 



(48) 



n— — oo 



Comparing ( [48| ) with (|43|), we see that the semiclassi- 
cal and uniform versions of (p\po) will be asymptotically 
equal if F(X,n) ~ J n (X/2)i n for large A. That is indeed 
true since for A ^> n 2 , 



F(X,n 




Table |l| shows that for A = 10 semiclassical and RM val- 
ues of a n differ by less than 0.02 up to n — 3. In the 
opposite limit, for A < 1, even ao evaluated semiclassi- 
cally is completely off. 
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VII. SCATTERING FROM A CORRUGATED 
WALL 

Let us apply the replacement-manifold method to 
a physical problem generating loop structure in phase 
space. Scattering of a plane wave by a corrugated wall 
has exactly such a property. This system has been used 
to model elastic scattering of atoms by solid surfaces (for 
review see JL2|| ). Consider a two-dimensional coordinate 
space divided into two parts by a periodically curved 
boundary consisting of the set of points with coordinates 
related by 







sin (3x 



(50) 



where (3 gives the spatial frequency of corrugation and e 
the maximum slope ^§ of the wall. A plane wave with 
momentum Kki is incident from the upper half-plane at 
an angle ai from the y axis, so the incident wave vector 
is ki = fc (sinai, — cosai). Classically, the wave reflects 
specularly from the curved boundary. The reflected rays 
are shown in Fig. pi Notice also the clearly visible caus- 




FIG. 8. Ray picture for a two-dimensional scattering from 
a corrugated wall. Only reflected rays are shown. 

tics. The Poincare surface of section in Fig. || displays de- 
pendence of momentum component p2x of reflected rays 
on coordinate i at a given distance yi from x axis. We 
can see familiar loops with constant area and predict the 
failure of semiclassical approximation when this area gets 
smaller than h. Below we uniformize the semiclassical 
solution using RMs and compare the outcome with the 
exact and primitive semiclassical results. 




FIG. 9. Original (solid line) and replacement manifolds 
(dashed lines) on the Poincare surface of section for scattering 
from a corrugated wall. 



There exist several approaches to solve this scatter- 
ing problem exactly. Garcia and Cabrera jl3| have thor- 
oughly compared merits of various methods. The main 
issue is the solution's convergence, which improves with 
decreasing corrugation parameter e. Luckily, this regime 
overlaps with the small-loop limit in which we are in- 
terested. We can therefore use the most straightforward 
approach based on the so-called Rayleigh hypothesis fLlf] 
which has been shown to be valid for e < 0.448 by Millar 
|l5| . Details of the solution are provided in the Appendix. 
A probability density plot of the wave function is shown 
in Fig. 0. 

The whole scattering problem can be formulated us- 
ing an analogue to the smooth-potential Lippmann- 
Schwinger equation. Invoking the Green's theorem, it 
can be shown that the total wave function satisfies 



dx \r)n(x) ■ V^ to t(r) (51) 



boundary 



where Go is the free-space Green's operator and h{x) is 
a normal unit vector at the boundary pointing into free 
space. Now we could use periodicity of the wall, expand 
V'tot(r) in its Fourier modes, 



Vw(r) = Vw(r) + ipscattr) 



2tt 



1 

2^ 



a « ' 



,ik2,' 



(52) 



and recover the linear system (A6) for coefficients a n . 
Our goal, however, is to make a connection with a 
semiclassical picture in the Poincare surface of section 
(Fig. [j]). That figure implicitly assumes that classi- 
cally only a single scattering event takes place before a 
ray leaves the wall permanently. If corrugation is deep 
enough (e > £ max , where e m ax < 1), multiple scatter- 
ing will occur for any incident angle a\. The larger e 
gets (within the range < e < e m ax), the smaller the 
maximum allowed incident angle a\. Let us consider a 




FIG. 10. Probability density plot for the exact quantum 
solution. Plot of the RM solution is indistinguishable to the 
eye. For detailed comparison, see Fig. [l2|. 

wall with very shallow corrugation (e = 0.2) and a small 
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angle of incidence (ai = 0.2), a situation which classi- 
cally allows only a single reflection. We also choose the 
incident-wave length A = ^ to be 0.3 times the period 
of corrugation. The surprising result we shall obtain be- 
low is that although there exist seven real Bragg peaks in 
the scattered wave, the semiclassics break down while the 
single-scattering approximation using replacement man- 
ifolds works with excellent accuracy. 

The simple semiclassical wave function may be evalu- 
ated e.g. by tracing of individual rays and employing the 
Van Vleck propagator |[(|. The action must be adjusted 
by correct Maslov indices, corresponding to reflection 
from a hard wall and to passage through caustics [ p7[ . 
A probability density plot of the semiclassical wave func- 
tion is shown in Fig. [n]. Note the caustics, which clearly 
separate regions with one, three, and five contributions 
to the scattered wave. Between caustics, solution looks 
qualitatively the same as the exact quantum analogue in 
Fig. By looking at a surface of section at y 2 = Air/ (3 
(where there are already at least three contributions for 
any x) in Fig. ^, we expect that semiclassical and quan- 
tum wave functions should disagree everywhere. This 
guess is confirmed in Fig. nJL 




FIG. 11. Probability density plot for the semiclassical so- 
lution. 

Now let us find a uniform expression for the scattered 
wave. The Lippmann-Schwinger equation ( |5l| ) is an in- 
tegral equation, and it is not immediately obvious how it 
could be used to find ipscat from ipinc ■ We circumvent this 
problem by proceeding in four steps: first we canonically 
transform the incident wave function to a new coordinate 
system r' in which the wall becomes straight. We solve 
the scattering problem in these coordinates, since there 
the semiclassical approach is exact. Then we transform 
back to the original coordinate system in which we prop- 
agate the wave using the free-space Green's function to 
obtain our final answer. Semiclassical approach is used 
for each separate step but their combination is evaluated 
exactly. 

The wall becomes flat if we apply a canonical transfor- 
mation 



x, y 



■ sin (3x 



generated by 



F 2 (r',k) = k x x' + k y {y' + ~ sm{3x'). 



(53) 



The incident wave function ipmcO^-i) can be transformed 
into new coordinates using 



<r'|k) sc = : le^( r '< k ) 

Z7T 



(54) 



In primed coordinates, the semiclassical solution of scat- 
tering is exact and "angle of incidence" equals "angle of 
reflection". Therefore h' ■ V?/>t t(r') may be replaced by 
2n' ■ Vtpi nc (r') in equation (pl|). The non-diagonal part 
of the semiclassical scattering matrix becomes 



f = 5-1 = 2 / dx'\r')n' ■ V'(r' 



(55) 



y'=0 



where 



y' . Transformation back to original coor- 



dinates is achieved by a complex conjugate of fl54|). The 
free-space propagation is accomplished by a semiclassical 
Green's function, exact in mixed representation 



(k 2x ,y 2 \G s c \k) = 



5(k 2x - k x ) e 



2ir(k 2 - k\ + ie) 



(56) 



Putting all four pieces together, we obtain an expression 
for the scattered wave on a surface of section, 



k 2 „y2\G s c T\ki) 

2 fd 2 k f dx' (k 2x , y2 \G^\k)(k\r') d y ,(r'\k) 

d 2 k- 



(57) 



2ik ly f ^ S(k 2x - k x ) e lk ^ 2 



(2tt) 5 / 2 



2ir(k 2 - kl + ie) 



x J dx' exp[i(ki x — k x )x' + i(k\ y — k y )— sin /3a;'] 

Before we proceed, it should be noted that the outcome 
of this four-step process, an integral representation of the 
scattered wave, can be viewed as a continuous superpo- 
sition of waves emanating from sources along the bound- 
ary, with strength proportional to the normal derivative 
of the incident wave. It can be shown, in fact, that this 
process is equivalent to a generalized Kirchhoff diffrac- 
tion method @. 

Returning to expression (^), we recognize that the 
exponent in integral over x' corresponds to a classical 
manifold 



k\ x + (ki y — k y )e cos fix 1 



(58) 



which is "begging" to be replaced by partial manifolds 
k x , n — k\ x + n[3x' with weights J„(A) because 



|A| = \k ly -ky\~ < 



'/3 



2ke 

T 



< i 







for the classically allowed momenta. The integral over x' 
is then simple to evaluate and is equal to 

p CO 

dx' Jn (X) exp[i(k lx + nP - k x ) x'} (59) 

n=— oo 
00 

= 2vr MX) S(k lx + nP - k x ) 

n=— 00 

Using this result and evaluating the trivial integral over 
k x , expression (|57| ) becomes 



k 2x ,y2\G s c T\k 1 ) 



00 . 

X! / dk y 



2ik 



' (2tt) 3 / 2 

e ik y V2 



n— — oo 



2n(k 2 - fe| + ie) 



(60) 



— ^n(A) 5(fc 2:c - k 2x . n ) 



where ki x . n = k\ x + nf3. Integral over k y picks up a pole 
at k y = k2y : n = \ k 2 — k\ x n and the final answer is 



(fc 2 „2/ 2 I Gq c T |ki 



E 



= fC2 B ,nJ/2 



X 



(61) 

k2x,n) 



v 2y ,n 



where A„ = — fc2y,n)§- The scattered wave is thus 
a superposition of traveling and evanescent waves with 
wavevectors k 2 n with coefficients 



7 — — «/n(An)- 



(62) 



These turn out to be very close to the exact quantum 
coefficients obtained in the Appendix, but here we have 
avoided the solution of a linear system (A6). 

The uniform expression (^2|) for coefficients of expan- 
sion ( |52] ) is surprisingly accurate. A two-dimensional 
density plot in Fig. [fojdoes not reveal any difference from 
the exact quantum solution while an analogous semiclas- 
sical plot clearly shows the caustics (see Fig. |ll|). Even 
if we look at the Poincare surface of section in a region 
with many classical caustics, the quantum and uniform 
solutions agree while semiclassical solution fails miser- 
ably (see Fig. |l2"| ). The uniform solution agrees with one 
found by Garibaldi et al. ]l9| and is somewhat more accu- 
rate than the same result without the kx y /k2y, n prefactor 
obtained by Hubbard and Miller |Io| by SCP approxima- 
tion and by Garibaldi et al. Jll| . As Garibaldi et al. (who 
obtain three different solutions differing only by prefac- 
tors) point out, these prefactors "should not be taken too 
seriously." We agree: for traveling modes, they do not 
cause large errors, and since all these solutions neglect 
multiple scattering, we cannot expect high accuracy of 
the already small coefficients of evanescent modes. 

As promised in Section IV, we have obtained almost 
identical result using RMs as Hubbard and Miller with 
SCP approximation [fill. They found their answer for a 



hard wall as a limit of the solution for a smooth potential 
due to Lennard- Jones and Devonshire Q . Although we 
were lucky in this case to obtain exactly solvable integrals 
for RMs, we can quite generally expect that the replace- 
ment manifolds will be much smoother than the original 
one, enabling their simple evaluation by SP method. 

We discuss the specific choice of parameters. Select- 
ing e < 1, a < 1 was necessary for the validity of the 
single-scattering approximation. At first, it appears that 
the simple semiclassical approximation should be accu- 
rate for this regime of very shallow corrugation, since 
scattering from a flat wall has an exact semiclassical an- 
swer. The correct criterion, however, is based on the area 
of loops in phase space (see Fig. ||). In this case, for small 
e and a, the height of loops is equal to the momentum 
kick in region of maximum slope, approximately ps and 
the width is half the period of corrugation, |j. The loop 
area is smaller than pe^ — irkdh (d is the depth of cor- 
rugation) and has to be compared to h. We arrive at 
a surprising result that the validity of semiclassical ap- 
proximation has nothing to do with the periodicity of 
corrugation, but only depends on the product kd of wave 
vector and corrugation depth. If this product is small, 
semiclassical approximation breaks down. In the oppo- 
site case of deep corrugation, semiclassical approximation 
works, although we might have to take into account clas- 
sical trajectories corresponding to multiple scattering. 




FIG. 12. Wave function on the surface of section at 
y = 4ir//3. Comparison of the exact quantum (dashed line), 
RM (solid line), and semiclassical (dotted line) wave func- 
tions. 



VIII. CONCLUSION 

We have successfully uniformized several intimidating 
situations which can arise when repeated small areas 
or loops arise from enclosing classical manifolds using 
the replacement-manifold idea. This notion can rejuve- 
nate physical intuition about the quantum wave func- 
tion and also be a very convenient approach to an accu- 
rate result. Two physically motivated cases were consid- 
ered (the homoclinic-like oscillation and the corrugated 
wall) along with several other contrived cases to test the 
method. The present method should work in more com- 
plicated time-dependent problems if used with Miller's 
IVR applied to Van Vleck's semiclassical propagator. Af- 
ter replacement manifolds are identified, the perturbative 
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expansion contributions to the wave function could be 
found by SP approximation applied to these manifolds. 

Future work along these lines includes the attempt to 
uniformizc a full homoclinic tangle. It remains to be seen 
how convenient and generally applicable the replacement 
manifold idea is. However, it seems clear that we now 
have a new method which simplifies many problems while 
reinstating an intuitive foundation. 
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Recalling the identity e ^m{(x) _ ^ = _ qo J m (\) e imx 
and denoting Ai = \ 2 . n — jjk 2yn , Q — 0x, the 

boundary condition may be written as 



E 



m— — oo 



Jm(Ai) e l 



^ 0.11 Jm(^2,r 



i(n+m)(> 



= 0. 



(A4) 



Replacing m by p = n + m in the second ter m and by 
p ee to in the first term in the brackets in (A4), we get 



E - 

p— — oo 



Jp(M) + J]] 0>nJp-n[M 



0. (A5) 



This equation will be satisfied for all x = % only if coef- 
ficients of e lp9 on the left-hand side vanish, 



APPENDIX A: EXACT SOLUTION OF 
SCATTERING FROM A HARD SINUSOIDAL 
WALL 

From conservation of energy, the scattered and incident 
waves must have the same magnitude of momentum. Due 
to the periodicity of the wall in the x direction, the re- 
component of momentum is constrained by the Bragg 
condition 

hx.n = hx + n/3, n e Z. (Al) 
Corresponding y-component of momentum k 2y ^ n = 
^Jk 2 — k 2x n will be real only if \ki x + nf3\ < k. Oth- 
erwise the reflected wave becomes evanescent: it decays 
exponentially with distance from a;-axis because k 2y is 
imaginary. To complete the solution, coefficients a n in 
the momentum Green's function 



J p (Xi 



^ a n Jp-n(A 2 , n ) = 0, VpeZ. (A6) 



We recognize (A6) as an infinite system of linear equa- 
tions for coefficients a n . In matrix form, A • a = b, with 

Aij ee Ji_j (j3\Jk 2 - (k lx + j/?) 2 )and h = -Ji{^k ly ). 

Although this system cannot be solved analytically, ex- 
pressing the solution in this form is a great improvement 
compared to what can be done for a general Dirichlct 
problem satisfying tim e- independent Schrodinger equa- 
tion. The system (A6) can be approximately evaluated 
numerically by restricting n and p to lie between limiting 
values ±ri max . 



(k2x\ki x ) = a n S(k 2x - k 2x . n ) S(k 2y - k 



must be found. This is accomplished by imposing 
the Dirichlet boundary condition in coordinate space. 
Namely, the total wave function iptot(r) = Vw(r) + 
fpscat(^) has to vanish along the corrugated wall. Using 



^ inc (r) = (r|k x ) - ^e 4kl ' r , 

Z7T 



(A2) 



f 1 °° 

Vw(r) = / d 2 k 2 (r|k 2 )(k 2 |ki) = — a n e 



definition of boundary (|50| ) and ( |Al|) for k 2x ^ n , the 
boundary condition becomes 



exp 



i{k\ x x + kiy— sin fix) 



^2 a„exp 



71— — OC 



i(ki x + n(3)x + ik2 y ,nTi sin (3x 



(A3) 
= 0. 
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